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We present an overview of recent developments in the detection of light bosonic 
dark matter, including axion, pseudoscalar axion-like and scalar dark matter, 
which form either a coherently oscillating classical field or topological defects 
(solitons). We emphasise new high-precision laboratory and astrophysical mea¬ 
surements, in which the sought effects are linear in the underlying interaction 
strength between dark matter and ordinary matter, in contrast to traditional de¬ 
tection schemes for dark matter, where the effects are quadratic or higher order in 
the underlying interaction parameters and are extremely small. New terrestrial 
experiments include measurements with atomic clocks, spectroscopy, atomic and 
solid-state magnetometry, torsion pendula, ultracold neutrons, and laser interfer¬ 
ometry. New astrophysical observations include pulsar timing, cosmic radiation 
lensing, Big Bang nucleosynthesis and cosmic microwave background measure¬ 
ments. We also discuss various recently proposed mechanisms for the induction 
of slow ‘drifts’, oscillating variations and transient-in-time variations of the funda¬ 
mental constants of Nature by dark matter, which offer a more natural means of 
producing a cosmological evolution of the fundamental constants compared with 
traditional dark energy-type theories, which invoke a (nearly) massless underly¬ 
ing field. Thus, measurements of variation of the fundamental constants gives us 
a new tool in dark matter searches. 


1. Introduction 

Dark matter (DM) remains one of the most important unsolved problems in con¬ 
temporary physics. Observations of stellar orbits about galactic centres from as 
early as the 1930s [1, 2], which were later rehned in the 1970s [3, 4], have indicated 
that the orbital velocities of stars remain approximately constant at large distances 
from the galactic centre (purple line in Fig. 1), rather than follow the Kepplerian 
dependence v oc If^/r (pink line in Fig. 1), which is expected from the observation 
that most stars are concentrated in the galactic core. These observations provide 
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strong evidence for the presence of DM in galaxies, which is predominantly located 
at moderately large distances away from the galactic centre. DM is a non-luminous, 
non-baryonic form of matter than interacts very weakly with itself and Standard 
Model (SM) matter. Observations of stellar orbital velocities in our local galactic 
neighbourhood give the cold (non-relativistic) DM energy density within our local 
galactic neighbourhood of [5]: 

pS = 0.4 GeV/cm^ (1) 

Further evidence for the existence of DM comes from gravitational lensing obser¬ 
vations of the Bullet Cluster [6-8], angular fluctuations in the cosmic microwave 
background (CMB) spectrum [9], and the need for non-baryonic matter to explain 
observed structure formation [10]. The latest Wilkinson Microwave Anisotropy 
Probe (WMAP) observations give a present-day mean DM energy density of [5]: 

PDM = 1.3 X 10"® GeV/cm^ (2) 



Fig. 1. Observed (purple line) and expected (pink line, which has the Kepplerian dependence 
V oc ^jy/r at large distances) orbital velocities of stars as a function of distance from the galactic 
centre. The discrepancy between the two radial functions is consistent with the presence of dark 
matter haloes in galaxies. 


In order to explain its observed abundance, it is reasonable to expect that DM 
has non-gravitational interactions with ordinary matter. Despite the overwhelming 
evidence for its existence, direct searches for DM via non-gravitational interactions 
with ordinary matter have not yet produced a strong positive result, leaving the 
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identity and non-gravitational interactions of DM in a state of mystery. Various the¬ 
oretically well-motivated candidates for DM have been proposed, including weakly 
interacting massive particles (WIMPs), axions, and weakly interacting slim parti¬ 
cles (WISPs). We refer the reader to the comprehensive reviews in Refs. [5, 11-13] 
for an overview of the theoretical motivation behind the main candidates, their role 
in cosmology and the main searches for these candidates. In traditional searches for 
WIMP DM (see e.g. Refs. [14-19]), which look for the scattering of WIMP DM off 
nuclei (Fig. 2), the sought effect is quartic in the underlying interaction parameters 
e' that parametrise the interaction between DM and nucleons: 

( 3 ) 

The smallness of the interaction parameters e' make further progress in these 
searches for WIMP DM very challenging. 



Fig. 2. Traditional WIMP dark matter detection experiments search for the scattering of WIMP 
dark matter off nuclei. The scattering cross-section associated with this process scales as: crgcat (x 




/M2)2 


In recent times, there has been a growing interest to use atomic and related sys¬ 
tems to directly search for DM. There is very strong motivation for the use of such 
systems, which to date have been employed with great success as high-precision fre¬ 
quency standards, in tests of the SM and as sensitive probes of new physics beyond 
the SM [20-24]. Atomic clocks are one of the most precise instruments every built 
by mankind, with the best current fractional inaccuracies of the order 10“^® [25-27]. 
Experiments with atomic Hg provide the most precise limits on the electric dipole 
moment (EDM) of the proton, quark chromo-EDM and P,T-odd nuclear forces, as 
well as the most precise limits on the neutron EDM and Quantum Chromodynam¬ 
ics (QCD) 6 term from atomic or molecular experiments [28] (ultracold neutron 
experiments give the best limits for the latter parameters [29]), while experiments 
with molecular ThO provide the most precise limit on the electron EDM [30]. Mea¬ 
surements and calculations of the 6s-7s parity nonconserving (PNC) amplitude in 
atomic Cs stand as the most precise atomic test of the SM electroweak theory to 
date, see e.g. Refs. [31-34], and are competitive with direct searches performed at 
hadron colliders [5, 35]. Experiments with atomic co-magnetometers [36-40], torsion 
pendula containing spin-polarised electrons [41, 42], and ultracold neutrons [43] pro- 
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vide some of the most stringent limits on CPT- and Lorentz-invariance-violating 
physics. Laser interferometer experiments have set the most stringent limits on 
gravitational wave detection to date [44-46]. 

In the present book review, we present an overview of recent developments in 
the detection of light bosonic dark matter, including axion, pseudoscalar axion-like 
and scalar dark matter, which form either a coherently oscillating classical field or 
topological defects (solitons), using a variety of high-precision laboratory measure¬ 
ments. We particularly emphasise new measurements, in which the sought effects 
are linear in the underlying interaction strength between dark matter and ordinary 
matter, and are easier to search for than traditional quartic effects such as those 
shown in Fig. 2. New terrestrial experiments include measurements with atomic 
clocks, spectroscopy, atomic and solid-state magnetometry, torsion pendula, ultra¬ 
cold neutrons, and laser interferometry. Astrophysical observations, such as Big 
Bang nucleosynthesis (BBN) and CMB measurements, assist in terrestrial searches 
by ruling out regions of the relevant parameter spaces. New astrophysical observa¬ 
tions that involve pulsar timing and cosmic radiation lensing can also be used to 
directly search for DM. We also discuss various recently proposed mechanisms for 
the induction of slow ‘drifts’, oscillating variations and transient-in-time variations 
of the fundamental constants of Nature by DM, which offer a more natural means 
of producing a cosmological evolution of the fundamental constants compared with 
traditional dark energy-type theories, which invoke a (nearly) massless underlying 
field. 

2. Axions 

2.1. The strong CP problem and the QCD axion 

When the SM was been developed during the 1970s, it quickly became apparent 
that there was an issue in the QCD sector as far the combined charge-parity (CP) 
symmetry was concerned. The QCD Lagrangian contains the P,CP-violating term 
[47-50]: 

'Cqcd = ( 4 ) 

where 9 is the angle that quantifies the amount of CP violation within the QCD 
sector, jpK = 14.5 is the colour coupling constant, and G and G are the gluonic 
field tensor and its dual, respectively. Account of weak interaction effects results 
in a shift of 9 from its bare value to the observable value 9 [51]. The angle 9 may 
in principle have assumed any value in the range —tt < 9 < -I-tt, but its observed 
value from measurements of the permanent static neutron EDM is constrained to 
be 10] < 10“^° [29]. The smallness of the observed value of 9 constitutes the strong 
CP problem. An elegant and the most widely accepted resolution of the strong CP 
problem was proposed by Peccei and Quinn [52, 53], in which the 9 parameter was 
interpretted as a dynamical field (the massive pseudoscalar axion, a): 9 ^ a{t)/fa, 
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where fa is the axion decay constant. Initially, the axion field is constant (0^1 
in the absence of fine-tuning of the vacuum misalignment angle 0i), but for times 
when nia ^ H, where H is the Hubble constant, the axion undergoes oscillations 
about the minimum of its potential (Fig. 3), which corresponds to d = 0, hence 
alleviating the strong CP problem [54-56]. 



Fig. 3. The QCD axion oscillates about the minimum of its potential, which corresponds to 0 = 0, 
thereby alleviating the strong CP problem. The frequency of oscillation is set by the mass of the 
QCD axion ma- 


2.2. Axions as cold dark matter 


Although the original PQWW model of the axion [52, 53, 57, 58] was quickly ruled 
out experimentally, the KSVZ [59, 60] and DFSZ [61, 62] models of the QCD axion 
turned out to be compatible with all terrrestrial and astrophysical observations 
(for some of the more recent invisible axion models based on the Peccei-Quinn 
symmetry, we refer the reader to Refs. [63-66]). The properties of the QCD axion 
are predominantly determined by the axion decay constant fa- In particular, the 
QCD axion mass ma is related to fa via the relation 

'10“ GeV 


m„ 


6 X10"^ eV 


fa 


( 5 ) 


For very weak couplings (i.e. for very large values of fa), axions are produced 
non-thermally in the early Universe. At temperatures well above the QCD phase 
transition, the QCD axion is effectively massless and the corresponding field can take 
any value, parameterised by 0i. The axion develops its non-zero mass ma (due to 
nonperturbative effects) when the temperature T < GeV, and for times when ma 3> 
H, the axion undergoes oscillations about the minimum of its potential at a = 0 
(Fig. 3). The resulting axion energy density, produced via this vacuum misalignment 
mechanism, is given (in terms of the critical energy density) by [54-56]: 

1.18 




0l 


fa 


1012 GeV 


( 6 ) 


For 01 ~ 1, axions saturate the present-day GDM content if fa ~ IQi^ GeV. For 
fa ^ 10i2 GeV, axion production via the vacuum misalignment mechanism would 
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have led to the overclosure of the Universe unless ^ 1, which may arise due to 
fine-tuning of the vacuum misalignment angle or anthropic selection [67-69]. We 
also note that the density of axions produced in the early Universe depends on the 
order of the cosmological events, in particular whether the Peccei-Quinn symmetry 
is broken prior to or following cosmic inflation. In the latter case, there may be 
additional contributions to the axion density of the same order as in Eq. (6) from 
the formation and decay of axionic topological defects, such as cosmic strings and 
domain walls [70]. 

Axions produced by the vacuum misalignment mechanism are very cold with 
almost no kinetic energy. Furthermore, if they are sufficiently light and weakly 
interacting, then these axions may survive until the present day and reside in the 
observed galactic DM haloes. If ma < 2me, then the axion lifetime is determined 
by its two-photon decay channel [12]: 


T(a —> 27) 


fl 


( 7 ) 


where jC'.yj ~ 1 is a model-dependent coefficient. For TOq ^ 24 eV, the axion lifetime 
exceeds the present age of the Universe. Thus, ultralight (sub-eV mass) axions, as 
well as axion-like pseudoscalar particles (ALPs) and scalar particles, for which no 
predictive mass formula akin to Eq. (5) exists, are good candidates for cold DM. 
Ultralight spin-0 bosons are good candidates for the dominant contributor to cold 
DM to very low particle masses. The simplest model-independent lower limit arises 
from the requirement that the de Broglie wavelength of the DM particles not exceed 
the halo size of the smallest galaxies, giving ma ^ 10“^^ eV. This simple estimate 
is in fact in good agreement with more rigorous limits obtained from cosmological 
and astrophysical investigations. Ultralight spin-0 DM would have inhibited cosmo¬ 
logical structure growth [71] in conflict with Lyman-alpha observations [72], unless 
^ 10“^^ eV. Ultralight spin-0 DM would have suppressed high-redshift galaxy 
formation, contrary to observations unless TOq > 10“^^ eV [73], while CMB obser¬ 
vations necessitate ma > 10“^^ eV [74]. We stress that these constraints only apply 
for ultralight spin-0 DM which is the dominant contributor to cold DM. Due to its 
effects on structure formation, ultralight spin-0 DM in the mass range 10“^^ —10“^° 
eV has been proposed [71, 75-77] to solve several long-standing astrophysical puz¬ 
zles, such as the cusp-core, missing satellite, and too-big-to-fail problems [78] (see 
also the earlier work of Ref. [79]). 

We note in passing that, while most interest resides in cold, non-relativistic, 
ultralight spin-0 bosons, the possibility of relativistic spin-0 bosons in the early 
Universe has also also been investigated. Relativistic particle species increase the 
total energy density in the early Universe (prei oc (1 -I- z)"^) and in turn increase 
the rate of cosmic expansion, which is parametrised by H (non-relativistic species, 
for which Pnon-rei OC (1 -f z)^ , do not affect the total energy density appreciably at 
early times). The presence of additional beyond-the-SM relativistic particle species, 
therefore, increases the neutron-to-proton ratio at the time of weak interaction 
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freeze-out, n/p = by causing freeze-out of the weak interactions, 

p + e~ ^ n + Ve and the corresponding crossing reactions, to occur at an earlier time 
(which corresponds to a larger freeze-out temperature Tp). Measurements and SM 
calculations of the primordial ^He abundance allow for an additional, relativistic, 
neutral spin-0 particle during BBN, while an additional relativistic spin-1/2 or spin- 
1 particle is excluded (the combination of the "^He abundance and the CMB value 
for the baryon-to-photon ratio p do not alter this conclusion) [80]. 

The number density of ultralight spin-0 fields per de Broglie volume readily 
exceeds unity, na/A^g 1. As a result, these bosons readily form a coherently 
oscillating classical field 

a(t) ~ ao cos(mat), (8) 

with an amplitude oq — V^Paxion/wa, where Paxion is the energy density associated 
with the bosonic field, and with a very well-defined oscillation frequency set by the 
boson mass. Over time, gravitational collapse of ultralight bosons into galaxies and 
their interaction with ordinary matter resulted in their virialisation, which led to a 
loss of perfect monochromaticity in their oscillation frequency 

(9) 

rUa 

where a virial velocity Uviriai ~ 10“^ is typical within our local galactic region. In 
the moving reference frame of our Solar System, a bosonic field has a non-zero 
average momentum and so the bosonic field takes the form 

a{r,t) = aQCOs{u}at-Pa-'r). (10) 


Despite the loss of perfect monochromaticity, the bosonic field remains coherent on 
time scales less than the coherence time 


”^coh 


277 



( 11 ) 


which is determined by the criterion that the additional phase accumulated over 
time in Eq. (10) due to virialisation remains less than 27r. 


2.3. Axion interactions and astrophysical constraints 


The axion couples to SM particles as follows (we consider only the couplings that 
are of direct interest to experimental searches): 




a 9 

fa 32^2 


GG 


G-yU 




( 12 ) 


where the first term represents the coupling of the axion field to the gluonic field 
tensor G and its dual G, the second term represents the coupling of the axion 
field to the electromagnetic field tensor F and its dual F, while the third term 
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represents the coupling of the derivative of the axion field to the fermion axial- 
vector currents and C/ are model-dependent coefficients. Typically, 

\C-y\ ~ 1 and \Cn\ ~ l^pl ~ 1 in models of the QCD axion [12, 81]. Within the DFSZ 
model, where the tree level coupling of the axion to the electron is non-vanishing, 
jCel ~ 1 [81]. However, within the KSVZ model, \Ce\ ~ 10“^, since the tree level 
conpling vanishes and the dominant effect arises at the 1-loop level [81]. For ALPs, 
the coefficients C-y and Cf are essentially free parameters, and the coupling to gluons 
is generally presumed absent. The common parameter in Eq. (12), which is of the 
most interest, is the axion decay constant fa- 

Astrophysical constraints on axion parameters greatly assist in laboratory 
searches for axions. For stellar axions, consideration of the axion production pro¬ 
cesses 7 -I- 7 —> 0 , 7 -|-e“—)'a-|-e“, 7 -|-A^—>'a-|-A^, the axion absorption processes 
a-|-e“ —>- 7 - 1 - 6 “, a + N ^ + N and the decay channel a —> 7 -I- 7 in stars gives 
the following lower bound on the mass of stellar axions [82]: 

^stellar ^ 25 keV, (13) 

while the requirement that stellar axion emission from helium-burning red giants 
not disrupt observed stellar evolution gives the stronger lower limit of [82, 83]: 

^stellar > 299 keV. (14) 

Energy loss from stars by solar axion emission (7 -I- 7 —> a) requires enhanced 
nuclear burning, which would lead to an increase in the solar ®B neutrino flux 
that contradicts observations unless the axion couples sufficiently weakly to the 
photon [84]: 


fa/C-y > 2 X 10® GeV. (15) 

The application of energy-loss arguments to the nucleon bremmstrahlung process 
N + N —>■ N + N + a in supernovae 1987A gives the following limit on the axion 
coupling to nucleons [84]: 


fa/CN > 10® GeV, (16) 

Energy loss mechanisms in stars through the Compton-like process 7 - 1 - 6 “ —> a -I- e“ 
and through the bremsstrahlung process e~ {Z, A) —> e“ -I- (Z, A) + a would 
excessively delay the onset of helium burning unless the axion couples sufficiently 
weakly to the electron [84]: 


la/Ca > 2 X 10® GeV, (17) 

with similar constraints from consideration of the increase in white dwarf cooling 
rates due to axion emission [84, 85]. 
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2.4. Traditional axion searches 

Haloscope and helioscope methods can be used to search for galactic and solar 
axions, respectively, via the axion’s coupling to the photon [86]. The traditional 
haloscope (ADMX) [87] and helioscope (CAST) [88] experiments (Fig. 4) have shed 
valuable light on our understanding on the possible axion parameter space for the 
axion-photon coupling. The lAXO helioscope experiment [89] will be the upgrade of 
the present CAST experiment. Searches for solar axions via the axio-electric effect 
with scintillator detectors have also been conducted [90-92]. Various ‘light-shining- 
through-wall’ [93-95] (Fig. 5) and vacuum birefringence [96, 97] searches for axions 
and ALPs via the axion-photon coupling have also been performed (for an overview 
of light-shining-through-wall searches for ALPs and various other light bosonic DM 
particles, we refer the reader to the review [98]). 



Fig. 4. Haloscope and helioscope experiments search for the conversion of galactic and solar 
axions, respectively, into photons in a strong applied magnetic field. The expected power generated 
by the conversion a ^ 7 scales as: Pa —»-7 c>c (l//a)^. 




Fig. 5. ‘Light-shining-through-wall’ experiments search for the transmission of photons through 
an impermeable material, due to the interconversion 7 ^ a ^ 7 in the presence of strong applied 
magnetic fields on either side of the barrier. The expected power generated by the interconversion 
7 ^ a ^ 7 scales as: P^^a^'y oc (l//a)^. 


Searches for ultra-light spin-0 bosons in tabletop experiments via the macro¬ 
scopic forces they would produce due to their couplings with the electron and 
nucleons (Fig. 6) have also been proposed [99]. Exchange of spin-0 bosons be- 
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tween fermions can produce either a spin-independent monopole-monopole poten¬ 
tial, a P, T-violating monopole-dipole potential with the cr • r correlation, or a fully 
spin-dependent dipole-dipole potential. Atomic magnetometry [100-107], torsion 
pendulum [108], differential force measurements [109] and ultracold neutron ex¬ 
periments [110-112] have collectively probed the axion-electron and axion-nucleon 
couplings over an expansive range of axion masses (see also Ref. [113] for con¬ 
straints on long-range interactions between spin-polarised geoelectrons deep within 
the Earth and the spin-polarised electrons and nucleons in laboratory experiments, 
and Ref. [114] for constraints on axion-electron and axion-nucleon interactions from 
a combination of terrestrial equivalence principle tests and astrophysical energy-loss 
bounds). Constraints on axion interactions through their mediation of spin-spin 
couplings in atomic systems have also been derived [115]. 



Fig. 6. Tabletop experiments can search for the effects of new macroscopic forces mediated by 
the exchange of ultralight spin-0 bosons. The induced energy shift due to such a new force scales 
as: 5e oc (l//a)^. Two vertices of scalar form ggcp'ipip produce a spin-independent monopole- 
monopole force. A vertex of scalar form, combined with a vertex of pseudoscalar form gpcp-iiji'y^ip, 
produce a P, T-violating monopole-dipole force. Two vertices of pseudoscalar form produce a fully 
spin-dependent dipole-dipole force. 


2.5. New axion searches 

A number of new proposals to search for axions have been put forward over the 
recent years. These proposals may be partitioned into two broad categories: 

(I) New-generation searches, where the underlying axion-induced effects are lin¬ 
ear in the combination a/fa {ao/fa — 4 x 10“^® for the QCD axion, which obeys 
Eq. (5) and which also saturates the local cold DM energy density in Eq. (1)), and 
are intrinsically much larger than the effects in traditional searches of Sec. (2.4) and 
those in (II) below. These new-generation searches are outlined in Secs. (2.5.1) — 
(2.5.4) that follow. 

(II) Searches, where the underlying effects are quadratic or higher order in 
the combination of axion parameters a/fa- These proposals are summarised in 
Sec. (2.5.5) below. 



September 4, 2015 0:11 


World Scientific Review Volume - 9.75in x 6.Sin 


Stadnik'Book'Chapter'Sepl page 11 


Manifestations of dark matter and variations of fundamental constants in atoms and astrophysical phenomenall 


2.5.1. ‘Axion wind’ effect 


The ^ = 1,2,3 components in the coupling of the derivative of the axion field to 
the fermion axial-vector currents in Eq. (12) lead to the following non-relativistic 
Hamiltonian [116-118]: 


ifeff(t) 


E 

f=e,p,n 


Cftto 


sin(mat) Va ' 




(18) 


which implies that a spin-polarised source of particles interacts with the axion 3- 
momentum, producing oscillating shifts in the energy of the spin-polarised source 
(which are linear in aojfa) at two characteristic frequencies: wi ~ nia and UJ 2 = 
27r/rsidereai, where Tgidereai = 23.93 hours is the sidereal day duration. This is the 
‘axion wind’ effect, which may be sought for using a variety of spin-polarised sources, 
for example, atomic co-magnetometers, torsion pendula and ultracold neutrons. 

Distortion of the axion field by the gravitational field of a massive body, such 
as the Sun or Earth, results in an additional axion-induced oscillating spin-gravity 
coupling (oscillating gravi-magnetic moments): H'^g(t) (x (C/Oo/Za) sin(mot) o-f -r, 
which is directed towards the centre of the gravitating body [117]. Eor couplings of 
the axion to nucleons inside the nucleus, isotopic dependence {Cn Cp) requires 
knowledge of the proton and neutron spin contributions in experiments that search 
for the ‘axion wind’ effect and oscillating gravi-magnetic moments. The proton and 
neutron spin contributions for nuclei of experimental interest have been calculated 
in Ref. [119]. 


2.5.2. Transient ‘axion wind’ effect 

Apart from the classical fields that ultralight axions and other spin-0 fields may form 
(see Sec. (2.2)), ultralight bosonic DM fields may also form topological defects, 
which arise from the stabilisation of the DM field under a suitable self-potential 
[120-126]. Topological defects, which make up a sub-dominant fraction of DM, are 
believed to function as seeds for structure formation [127]. Eor some of the more 
recent developments on topological defects, we refer the reader to Refs. [128-132], 
while for the classical review, we refer the reader to Ref. [133]. 

While stable domain wall structures that consist of the QCD axion would lead 
to disastrous consequences in cosmology by storing too much energy [126], domain 
walls and other topological defect structures consisting of ALPs or scalars are viable 
for certain combinations of parameters. Topological defects that consist of ALPs 
may interact with fermion axial-vector currents via the /i = 1,2,3 components of 
the derivative coupling in Eq. (12), which in the non-relativistic limit reads [134]: 

f=e,p,n “ 

Eq. (19) implies that a spin-polarised source of particles will temporarily interact 
with a topological defect as the defect passes through the system. A global network 
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of detectors, such as atomic co-magnetometers, has been proposed to detect these 
correlated transient-in-time signals, produced by the passage of a defect through 
Earth [134, 135]. 


2.5.3. Oscillating P,T-odd electromagnetic moments 

Interaction of the QCD axion field with the gluon fields in Eq. (12) produces an 
oscillating neutron EDM [117, 136]: 

dn{t) ~ 1.2 X 10“^® ^ cos(mat) e • cm, (20) 

Ja 

which induces oscillating nuclear Schiff moments [117, 136] and oscillating nuclear 
magnetic quadrupole moments [137]. In nuclei, a second and more efficient mech¬ 
anism exists for the induction of oscillating electromagnetic moments by axions — 
namely, the P, T-violating nucleon-nucleon interaction that is mediated by pion ex¬ 
change, with the axion field supplying the oscillating source of P and T violation 
at one of the ttNN vertices [117] (Fig. 7). 




Fig. 7. Main process responsible for the induction of oscillating P,T-odd nuclear electromagnetic 
moments by an oscillating axion field. The black vertex on the left is due to the usual strong 
P,T-conserving ttNN coupling (Qt^nn = 13.5), while the magenta vertex on the right is due to 
the axion-induced P,T-violating nNN coupling {g^NN — 0.027ao//a. cos(maf)) [117]. 


Axion-induced oscillating P^T-odd nuclear electromagnetic moments can in turn 
induce oscillating EDMs in atoms and molecules. In diamagnetic species (J = 0), 
only oscillating nuclear Schiff moments (which require I > 1/2) produce an oscillat¬ 
ing atomic/molecular EDM (oscillating nuclear EDMs are effectively screened for 
typical axion masses, as a consequence of Schiff’s theorem [138]). Two atoms that 
are of particular experimental interest are ^®®Hg and ^^®Ra, for which the axion 
induces the following oscillating EDMs [117]: 

d(^®®Hg) = —1.8 X 10“^® ^ cos(mat) e • cm, 

Ja 

fi(^^^Ra) = 9.3 X 10“^^ ^ cos(mat) e • cm, 

Ja 

with the large enhancement in ^^^Ra compared with ^®®Hg due to both collective 
effects and small energy separation between members of the relevant parity dou¬ 
blet, which occurs in nuclei with octupolar deformation and results in a significant 


( 21 ) 

( 22 ) 
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enhancement of the nuclear Schiff moment [139, 140]. A possible platform to search 
for the oscillating EDMs of diamagnetic atoms in ferroelectric solid-state media has 
been proposed in Ref. [141]. 

Paramagnetic species (</ > 1/2) offer more rich possibilities. Firstly, axion- 
induced oscillating nuclear magnetic quadrupole moments (which require / > 1) also 
produce an oscillating atomic/molecular EDM [137], which is typically larger than 
that due to an oscillating nuclear Schiff moment (since magnetic quadrupole mo¬ 
ments are not subject to screening of the applied electric field by atomic/molecular 
electrons). Secondly, an entirely different mechanism exists for the induction of 
oscillating EDMs in paramagnetic species, through the derivative interaction of the 
axion field with atomic/molecular electrons in Eq. (12). The ^ = 0 component of 
this term mixes atomic/molecular states of opposite parity (with both imaginary 
and real coefficients of admixture), generating the following oscillating atomic EDM 
(due to the real coefficients of admixture) in the non-relativistic approximation for 
an S'i /2 state [117]: 


da{t) ~ 


Ceaomlas 
faa ® 


sin(mat), 


(23) 


where as is the static scalar polarisability. Fully relativistic Hartree-Fock atomic 
calculations are in excellent agreement with the scaling da o: as in Eq. (23) [137, 
142]. The imaginary coefficients of admixture in the perturbed atomic wavefunction 
produce P-violating, T-conserving effects in atoms, while the analogous imaginary 
coefficients of admixture in the perturbed nuclear wavefunction (due to the axion- 
nucleon interaction via the fx = 0 component of the third term in Eq. (12)) produce 
P-violating, T-conserving nuclear anapole moments [117, 137, 142]. 

An axion or ALP field may also induce oscillating EDMs in paramagnetic species 
via perturbation of the electron-nucleon Coulomb interaction by the axion-photon 
interaction of Eq. (12) [143] (Fig. 8). 



Fig. 8. Induction of an oscillating electric dipole moment in a paramagnetic atom or molecule 
may occur as a result of the perturbation of the Coulomb interaction of atomic/molecular electrons 
and nucleons by the axion electromagnetic anomaly. 
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2.5.4. Oscillating electric current flows along magnetic field lines 

Interaction of an axion field with the photon field in Eq. (12) causes an oscillating 
electric current to flow along magnetic field lines [144]. Conversely, in the presence 
of an externally applied magnetic field i?o, an axion field induces an electric current 
density 

3a = Bo sm{mat), (24) 

which in turn produces a magnetic field Ba that satisfies V x B^ = j and can 
be amplified with an LC circuit and then detected using a magnetometer [144]. An 
analogous strategy has also been proposed for the detection of hidden photons [145]. 

2.5.5. Other proposals 

Ref. [146] proposes to modify existing microwave cavity searches for axions through 
the insertion of radio-frequency cavities into dipole magnets from particle acceler¬ 
ators, wiggler magnets developed for accelerator-based advanced light sources, and 
toroidal magnets similar to those used in particle physics detectors, while Ref. [147] 
proposes to modify existing microwave cavity detectors using an open Fabry-Perot 
resonator. Ref. [148] suggests to search for the electromagnetic radiation emit¬ 
ted by conducting surfaces when they are excited by ALPs (or hidden photons). 
Refs. [149, 150] propose to search for a Shapiro step-like signal induced by axions 
in Josephson junctions. Ref. [151] proposes to exploit nuclear magnetic resonance 
to search for axion-mediated forces, while Ref. [152] proposes to search for axion- 
induced atomic transitions in macroscopic samples. We also note that an external 
time-dependent magnetic field may enhance the local energy density stored in an 
ALP field by several orders of magnitude [153], which may have application to 
laboratory axion searches. 

3. Variations of the fundamental constants of Nature 
3.1. Traditional observations and models 

The idea that the fundamental constants of Nature might vary with time can be 
traced as far back as the large numbers hypothesis of Dirac, who hypothesised that 
the gravitational constant G might be proportional to the reciprocal of the age 
of the universe, G (x 1/t (which was later shown to be inconsistent with observa¬ 
tions) [154-156]. Since Dirac’s initial hypothesis, a number of models, in which 
the fundamental constants vary with space and time, have been proposed and in¬ 
vestigated, including Bekenstein models [157-160], string dilaton models [161-163], 
chameleon models [164], and via quantum effects induced by cosmological renor¬ 
malisation group flow [165-168] (see also the review [169] for an overview of other 
models). Anthropic arguments point out that ‘fine-tuning’ of the fundamental con¬ 
stants is required for life to exist — if the fundamental constants were even slightly 
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different in our region of the Universe, then life could not have appeared. Variation 
of the fundamental constants of Nature across space provide a natural explanation 
to such ‘fine-tuning’: mankind simply appeared in a region of the Universe where 
the values of the fundamental constants are suitable for our existence. 

Many laboratory, terrestrial and astrophysical searches for possible variations in 
the fundamental constants have been conducted. In any search for variations of the 
fundamental constants, the observable must be dimensionless (for otherwise, the 
observable would depend on the choice of units). Atomic transition frequencies are 
sensitive to variations in the electromagnetic hne-structure constant, a = 
and particle masses; for instance, the leading-order relativistic corrections to the 
hydrogenic Dirac energy levels scale as oc (Za)^. Atomic clock and spectroscopy 
measurements in the laboratory using a wide range of systems have provided some 
of the most precise limits on temporal drifts in a and particle masses to date [170- 
185]. The most stringent laboratory limits on temporal variations of a come from 
Al+/Hg+ [182] and Yb^/Yb’*' [185] clock comparison experiments: 

\{da/dt)/a\ < 2 x 10“^^ year“^, (25) 

while the most stringent laboratory limits on temporal variations of the electron-to- 
proton mass ratio TOe/wp come from Yb“'"/Cs clock comparison experiments [184]: 

l(i(me/mp)/(me/mp)l < 2 x 10“^® year“^. (26) 

A variety of systems have been proposed to provide improved laboratory limits on 
temporal variations in the fundamental constants. Optical transitions in highly- 
charged ionic species that are near the crossing of different electronic configura¬ 
tions are very sensitive to variations in a [186, 187]. Molecules, in which there 
is near cancellation between hyperfine and rotational intervals [188], ground-state 
fine-structure and vibrational intervals [189], and omega-type doubling and hyper¬ 
fine intervals [190], have enhanced sensitivity to a, rue/mp and the ratio of light 
quark masses to the QCD scale w^/Aqcd- Molecular ions, in which transitions oc¬ 
cur between nearly degenerate levels of different nature, have enhanced sensitivity 
to a and melrrip [191]. Transitions in some nuclei, such as the ultraviolet transi¬ 
tion between the ground and first excited states in the ^^^Th nucleus, have highly 
enhanced sensitivity to variations in a [192]. Most recently, laser interferometers 
have been proposed to search for variations of a [193]. The phase accumulated in 
an interferometer arm, <1> = wLjc^ changes if the fundamental constants change 
(the arm length L = Nas, where ob is the Bohr radius, and atomic frequency w 
both depend on the fundamental constants), according to 5<1> ~ $ 5aja, with a 
typical accumulated phase of $ ~ 10^^ in a single passage for an optical transition 
and L = 4 km. Multiple reflections enhance the accumulated phase and effects of 
variation of a by the factor Neff ~ 100. 

Stringent terrestrial limits on the temporal variation of a have been determined 
from the Oklo phenomenon. Roughly 1.8 billion years ago, a uranium-rich natural 
nuclear reactor in Oklo went critical, consumed a portion of its fuel and then shut 
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down several million years later. The isotopic ratio (neither of which 

are fission products) in Oklo ores is much lower than in ordinary samarium, which 
can be interpretted as been due to the depletion of ^^®Sm by the capture of thermal 
neutrons by n+ ^^®Sm —)■ ^^^Sm + y, which is dominated by the capture of 

a neutron with an energy of about 0.1 eV. This low-energy resonance arises due to 
the near cancellation between the electromagnetic and strong interactions, and the 
position of this resonance depends strongly on a, from which the following limit on 
temporal variations of a is obtained [194]: 

\{da/dt)/a\ < 10“^^ year“^. (27) 

Shifts in quasar absorption spectral lines provide a powerful tool to search for 
spatial variations in a and mejmp [195-202]. The most recent independent data 
samples from the VLT and Keck telescopes both indicate the presence of a spatial 
gradient in a [201]: 

— ~ 10-1® ly-i. (28) 

a 

A consequence of this astronomical result is that, since the solar system is moving 
along this spatial gradient, there should exist a corresponding temporal shift in a in 
Earth’s frame of reference at the level 5a/a ^ 10“i® year”! [203]. Finding this vari¬ 
ation with laboratory experiments may independently corroborate the astronomical 
result. The dynamics of electron-proton recombination is governed by a and me- 
CMB measurements thus provide a means of probing possible variations of the fun¬ 
damental constants, with sensitivities at the fractional level of ~ 10“^ — 10“® for 
variations of a and me from the present-day values, respectively [204-209]. The pri¬ 
mordial light elemental abundances, which are produced during BBN, are sensitive 
to changes in the fundamental constants, with sensitivities to temporal variations in 
the constants from the present-day values at the fractional level ~ 10“^ [210-217]. 

The underlying cause of any possible variations in the fundamental constants 
is still an open question. Traditional dark energy-based models, which predict a 
cosmological evolution of the fundamental constants, invoke a (nearly) massless 
underlying field. DM-based models offer a more natural approach to producing 
variations of the fundamental constants of Nature, since the underlying DM fields 
do not necessarily need to be exceedingly light. In the remaining sections, we present 
an overview of the mechanisms through which slow ‘drifts’, oscillating variations and 
transient-in-time variations of the fundamental constants may be induced by DM. 


3.2. Scalar interactions and constraints 

Scalar fields may interact linearly with SM matter as follows 


A. 




/ 


y 


^lin. _ 

-^scalar 


4 


( 29 ) 
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where the first term represents the coupling of the scalar field to the electromagnetic 
field tensor F, the second term represents the coupling of the scalar field to the 
fermion bilinears //, while the third term represents the coupling of the scalar 
field to the massive vector boson wavefunctions. The various Ax that appear in 
Eq. (29) are the respective new physics energy scales (analogous to the axion decay 
constant fa in Eq. (12)), which are independent of one another and are known to 
be very large energy scales from fifth-force searches (the exchange of an ultralight 
scalar particle between two SM fermions produces an attractive Yukawa potential, 
V{r) = —m'jj4:TTA'j ■ e“'"'^’’/r), which include lunar laser ranging [218, 219] and the 
EotWash experiment [220, 221]. The most stringent bounds for the scalar masses 
^ lO”^"* eV are from the EotWash experiment [220, 221]: 

> 3 X 10^2 Qgv, Ae > 2 X 10^1 GeV, 

Aq = {rud + mu)AdAu/(mdAu + m„Ad) > 5 x 10^^ GeV. (30) 

Stellar energy-loss arguments applied to the bremmstrahlung processes e~ + "‘He —)■ 
e~ + ^He -I- (j) give the following limit on the linear coupling of a scalar field to the 
electron [ 222 ]: 

Ae > 4 X 10‘° GeV (31) 

while similar arguments applied to the Gompton-like process 7 -I- ^He (j)+ "‘He 
give the following limit on the linear coupling of a scalar field to the proton [ 222 ]: 

Ap > 2 X 10‘° GeV. (32) 


Likewise, (pseudo)scalar DM may also interact quadratically with SM matter as 
follows 


^quad. 

■^scalar 


{AfY 4 


E 



^///+E 


M 2 

iKY 2 


VaV\ 


(33) 


where the A'^ are much less severely constrained than the corresponding Ax, from 
astrophysical observations and fifth-force searches. Stellar energy-loss arguments 
applied to the photon pair annihilation process 7 -I -7 —t 4>+4> constrain the quadratic 
coupling of a scalar field to the photon [223]: 


A'.y > 3 X 10^ GeV, 


(34) 


while similar arguments applied to the nucleon bremmstrahlung process N + N ^ 
N + N + (j) + (j> give the following limit on the quadratic coupling of a scalar field 
to the proton [223]: 

A], > 15 X 10^ GeV, (35) 

and similarly for the bremmstrahlung process e~ + {A, Z) ^ e~ + (A, Z) + (j) + (f, 
which constraints the quadratic coupling of a scalar field to the electron as follows 
[223]: 


A'g > 3 X 10^ GeV. 


(36) 
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For the quadratic couplings of Eq. (33), a fifth-force is produced in the leading 
order by the exchange of a pair of ^-quanta, which generates a less efficient V (r) ~ 
—TOj/647r^(A)-)^ • 1/r^ attractive potential, instead of the usual Yukawa potential in 
the case of linear couplings. The resulting constraints from fifth-force experiments 
are hence weakened significantly [223]: 

A'p>2x 10^ GeV, (37) 

which hold for < lO”"* eV [224]. 


3.3. Oscillating variations of the fundamental constants 


The couplings in Eq. (29) alter the electromagnetic fine-structure constant a and 
particle masses as follows 


1 - (j)/A^ 


1 + f- 


5m f (p 6 My (p 

mf Af^ My Ay’ 


(38) 


and, likewise, the couplings in Eq. (33) alter a and the particle masses as follows 


a —)• 


a 

1 - <py{Ap)^ - 



5m f 

mf 


(p"^ 5 My 

WfF^ 



(39) 


The couplings of an oscillating scalar field to the SM fields via the linear interactions 
in Eq. (29) [193, 225, 226] and via the quadratic interactions in Eq. (33) [193, 
226, 227] produce oscillating variations of the fundamental constants, which can be 
sought for with high-precision terrestrial experiments involving atomic clocks [225- 
227] and laser interferometers [193]. A multitude of atomic, highly-charged ionic, 
molecular and nuclear systems can be used in clock-based searches, see the reviews 
[228, 229] for summaries of the possible systems. The first laboratory clock-based 
search for oscillating variations of a has very recently been completed [230], and 
the results have been used to place stringent constraints on the photon interaction 
parameters A.^ [230] (Fig. 9) and A'^ [226, 227] (Fig. 10) for the range of scalar DM 
masses 10“^"* eV < < 10“^® eV. 


3.4. ‘Slow’ drifts of the fundamental constants 


The coupling of an oscillating scalar field to the SM fields via the quadratic couplings 
in Eq. (33) produces not only oscillating variations of the fundamental constants, 
but also ‘slow’ drifts of the fundamental constants through the ((/)^) = (Pq/2 terms in 
Eq. (39) [226, 227], which is related to the local ambient density of non-relativistic 
scalar DM via the relation Pscaiar — The dynamics of electron-proton re¬ 

combination is governed by a and me. CMB measurements constrain the quadratic 
interactions of p with the photon and electron as follows [227]: 


A' > 


1 eV^ 

mf, 


(40) 
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Fig. 9. Region of scalar dark matter parameter space ruled out for the linear interaction of 
(p with the photon. Region below red line corresponds to constraints from atomic dysprosium 
spectroscopy measurements [230]. Region below black line corresponds to constraints from fifth- 
force experimental searches [220, 221]. 


Changes in the fundamental constants during and prior to BBN alter the primordial 
abundances of the light elements. The observed and calculated (within the SM) ratio 
of the neutron-proton mass difference to freeze-out temperature at the time of weak 
interaction freeze-out (tp ~ 1-1 s), which determines the abundance of neutrons 
available for BBN (with the vast majority of these neutrons ultimately being locked 
up in ^He), constrain the quadratic interactions of (f) with the photon, light quarks 
and massive vector bosons as follows [226, 227]: 


1 J 

f 0.08 1.59 

rud ruu 

3.32 4.65 ] 

ml] 

. (Aiy)^ ' md-rriu 

mr iK)^\ 

+ (AV)2 (Ay 2 ; 


~ (1.0 ± 2.5) X 10-20 eV-'‘, (41) 

when the scalar field is oscillating at tp « 1.1 s (m^ 10“^® eV). When the 

scalar field had not yet begun to oscillate at tp « 1.1 s (m^ ^ 10“^® eV), the 
corresponding constraints on the quadratic interactions of 4> are [226, 227]: 

1 / 7710 r 0.08 1.59 / TTid 777„ \ 3.32 4.65 ' 

ml\3x 10-10 eVy [(Ag2 rrid-mu iK)'^ ) ^ . 

~ (0.5 ± 1.3) X 10-20 eV-'‘, (42) 
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while the constraints on the linear interactions of (j) are [226]: 

1 / y/'^rO.OS 1.59 (rrid 3.32 4.65' 

V3 X 10-16 eVy _ md-mu\f^d A„ / Kw A^ 

~ (0.4 ± 1.0) X 10"“ eV^ (43) 

The constraints on the photon interaction parameter A!) from CMB and BBN mea¬ 
surements are shown in Fig. 10. 



-20 -15 -10 -5 0 

■»€) 

Fig. 10. Region of scalar dark matter parameter space ruled out for the quadratic interaction 
of (j) with the photon. Region below red line corresponds to constraints from atomic dysprosium 
spectroscopy measurements [226, 227]. Region below yellow line corresponds to constraints from 
CMB measurements [227]. Region below blue line corresponds to constraints from comparison 
of measurements and SM calculations of the ratio QnpITp [226, 227]. Region below black line 
corresponds to constraints from stellar energy loss bounds and fifth-force experimental searches 
[223]. 


3.5. Transient-in-time variations of the fundamental constants 

The coupling of a scalar DM field that comprises a topological defect with SM fields 
via either of the couplings in Eqs. (29) or (33) alters the fundamental constants and 
particle masses inside the defect, giving rise to local transient-in-time variations as 
a defect temporarily passes through this region [231, 232]. These transient-in-time 
variations can be sought for using a global network of detectors, including atomic 
clocks [231] and laser interferometers [193], as well as a network of pulsars [232], 
such as the international pulsar timing array [233], or binary pulsar systems [234]. 
For sufficiently non-adiabatic passage of a defect (a relatively thin and/or rapidly 
travelling defect) through a pulsar, a topological defect may trigger a pulsar glitch 
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[232, 234], Numerous pulsar glitches have already been observed (see e.g. Ref. [235] 
for an overview), but their underlying cause is still debated (see e.g. Ref. [236] for 
a review). 

4. Outlook 

Exciting developments in dark matter searches are expected over the next few years. 

Effects that are linear in the interaction constant between dark matter and ordinary 
matter provide strong motivation for a new generation of searches for ultralight ax- 
ion and scalar dark matter. The first such laboratory search for ultralight scalar dark 
matter by means of atomic spectroscopy measurements in dysprosium has already 
been completed, placing new constraints on the linear and quadratic interactions 
of ultralight scalar dark matter with the photon that surpass previous constraints 
by many orders of magnitude. A number of other laboratory searches for ultralight 
dark matter using atomic and solid-state magnetometry, atomic clocks, interferome¬ 
try, torsion pendula and ultracold neutrons are either already in progress or planned 
to commence in the near future. These experiments are expected to yield limits on 
the interaction parameters of ultralight dark matter with ordinary matter that are 
many orders of magnitude better than existing limits and, more importantly, offer 
reinvigorated hope for the unambiguous direct detection of dark matter. 
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